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^ _ Abstract 

I Graphene is characterized by chiral electronic excitations. As such it provides a perfect testing 

■ ground for the production of Klein pairs (electron/lioles). If confirmed, the standard results 
for barrier phenomena must be reconsidered with, as a byproduct, the accumulation within the 
barrier of holes. 

Oi 

I. INTRODUCTION 

^ , Graphene is a subject of intense interest in the scientific community. After the assignment of the 2010 

Nobel prize in Physics, for the experimental realization of graphene obtained by extracting single 
atomic thick crystallites from bulk graphite, to Geim and Novoselov [HIS], many papers appear every 
day with new and interesting discussions about the properties of this intriguing material. Due to this 
vast amount of literature periodic updates are needed to cover the recent progress. The most explored 
aspects of graphene physics are related to its electronic properties. Due to the fact that electrons 
, propagating through graphene lose their effective mass, they are described by a Dirac-like equation 

■ rather than the Schrodinger equation. This open the doors to new possibilities to explore effects 
that were previously inaccessible. The massless Dirac equation in graphite intercalation compounds 

psj ■ was first pointed out in 1984 [3]. Nevertheless, the first observation of the Dirac fermion nature of 

, electrons in grephene [lllS] was only provided when graphene was mechanically exfoliated on SiOs 

■ [lj[2] . Relativistic phenomena related to Zitterbewegung [6] , chirality and minimal conductivity [7l|8] , 
Klein tunneling [5] , quantum Hall effect [TOIIII] , Schwinger production [T^] , Casimir interactions [13] , 
transformations of discrete Landau level spectrum to a continuum of extended states in presence of 

, a static electric field [M] have been investigated and when possible tested by using this reach and 

■ surprising material. For an extensive theoretical review of the properties of electron transport in 
\ graphene, we refer the reader to the report of Castro [15] . 

In the last year, as observed by Geim |16] . the initial skepticism with respect to graphene applica- 
tions was gradually evaporating and graphene has rapidly changed its status from newcomer to star. 
The study presented in this paper was inspired by the possibility to understand, by an experiment 
with graphene, the physics which governs a well known phenomenon of the Dirac equation known 
as Klein paradox [17], i.e. for the step potential the reflected flux can be larger than the incident 
one [18]. This relativistic effect is attributed to the fact that sufficiently strong potentials open the 
door to pair production and consequently to positrons inside the potential region. This result is gen- 
erally accepted for a step potential but not always quoted for a barrier potential where plane wave 
solutions seem not to exhibit pair production jl9| . Pair production at the discontinuity of square po- 
tentials has never been observed experimentally because its observation requires very strong electric 
fields [20]. The transmission of electron through graphene heterojunctions [2TH23] could now provide 
as an experimental signature of the Klein pair production. 

Graphene is a two dimensional structure characterized by chiral (mass zero) electronic exitations 
one of whose equation of motion reads 

- ihvpO- -V tpi.ir) ^ Eippir) , (1) 
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where <t = ((Ti,cr2,0) and Vp « 10 m/sec is the Fermi velocity. This equation mimies one of the 2- 
dimensional free Weyl equations [T51[^ for a (Weyl) spinor if we substitute Vp by the velocity of light 
c. In the following, we shall present the results for the Weyl equation for a potential barrier via the well 
tried step method pSl - ET] . Consequently, we "reconsider" the hypothesis of Klein pair production [T7] 
for the step potential. The step potential is actually a purely abstract construct, because of its infinite 
extension (same abstraction of plane waves). However, as we have demonstrated elsewhere with the 
Dirac equation the step potential is a useful tool for general piece- wise potential calculations. 
Contrary to a common misconception, it is compatible with a very large barrier potential. We will 
have the opportunity to reassert this fact here with our results for the Weyl equation. Furthermore, 
and this is essential to our conclusions, if Klein pair production occurs for a step, it must also occur 
for a barrier potential |19| . 

For the Dirac equation [TH1I2H], the so-called Klein energy zone (where Klein pair production is 
possible) requires a step potential height Vq of at least 2 m {m the mass of the fermion). For a free 
electron this means Vg > 1 MeV. Not a typical laboratory potential energy step. Not surprisingly, 
experimental verification of Klein pair production has not been achieved to date. Neutrinos are 
a possible alternative choice of particle since they are now known to have a very small mass, but 
the weak nature of their interactions introduces formidable experimental difficulties |29] . The best 
alternative "particle" is that of electrons within a dielectric with effective masses much smaller 
than the free electron mass. Graphene, with a null effective electron mass, is the most promising 
material of all to study. A step potential of any height could test the creation of Klein pairs (electron- 
hole) and as we have anticipated this has consequences even for the barrier potential seen within the 
graphene structure. 

In the next section, we calculate, in 3-dimensions for completeness, the chiral interaction with a 
step and subsequently in section III that with a barrier potential. The arguments require a discussion 
of an incoming wave packet and resolves the apparent paradox that a step has a single transmitted 
wave packet (both for the diffusion zone and the Klein zone) while a long barrier gives rise to infi- 
nite wave packets. We shall again argue that nevertheless the step and barrier results are perfectly 
compatible. In section IV, we discuss Klein pair production. We also show that for this case the 
transmission/refiection amplitude series are divergent. The finite standard (matrix) results are radi- 
cally different. They do not involve pair production. The application to graphene is straightforward. 
We argue that the barriers within graphene act as a "well-trap" for the created "antiparticles" . This 
should result in a growth in time of positive charge within the well/barrier. We suggest the detection 
of such an accumulating charge as proof of Klein pair production. 



II. WEYL INTERACTION WITH A STEP POTENTIAL 

The spinor Weyl equation for a positive helicity (massless) state is [18] 

- zcr • V*(r,t) = ift *(r,t)/c . (2) 

Plane wave solutions of the kind tpir) = u{p) exp[i (p • r — Ei)/li\, where E = \p\c and p is the 
particle three-momentum, result in the spinor equation 

a ■ p u{p) ^ Eu{p) / c , (3) 

the solution of which is 

where the upper index t represents the transpose. Our step potential is chosen in the x variable. 
Without loss of generality, we set the discontinuity at x = 0, 

V{x) = { for a; < region I, VJ, > for x > region II } . 

In region I (a; < 0), we have set the incoming momentum as p = (^1,^2,^3) with 



E + p;c 



(4) 



In region II (x > 0), only pi changes. We call the new momentum q = (gijPajPs) with 



qic 



The continuity equation for 5'(r,t) at a; = is 
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After a straightforward calculation, we find 
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Starting from the reflection probability, 
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there are specific cases to be considered. 

2 

The situation when < 0, i.e. when pi + pi > [{E — Vo)/c] , implies an imaginary and 
whence to "tunneling" when penetration into the classically forbidden region II occurs. An imaginary 
momentum also implies Re[a] = 0, consequently |ro| — I. Wc arc not interested in evanescent 
solutions in this paper and henceforth consider only real momentum qi . 

Diffusion occurs when E > Vo + y^pl + pl c. In this case Re [a] > and consequently |ro| < 1. We 

c, we have again a real momentum qi but now 

This Klein 



\/vl + Vl 

treat this case in the next section. For E <Va 
E 



+ 

K +P3C < E—Vo + \Jp\ +pIc < and consequently Re [a] < which implies |ro| > 1. 
energy zone is thus characterized by an oscillatory (free) plane wave in region II. The corresponding 
flux could, in principle, be interpreted as an additional incoming particle beam from x = +00 hence 
resulting in an excess of "reflected" particles in region I, or in alternative as an antiparticle flux flowing 
in the positive x direction and "created" at a; = 0. The former choice contradicts our initial conditions 
of a sole incoming particle wave from the left (x = —00) it also implies the unpalatable concept of free 
particles that exist and propagate freely in the classically forbidden region below potential. The latter 
solution implies the creation of Klein pairs. The "antiparticles" (energy —E) travel to the right while 
the particles created (energy E) add to the reflected incoming particles and travel to the left, thus 
yielding |ro| > 1. This latter interpretation is consistent with the fact that antiparticles see a potential 
of —Va and hence are above potential free antiparticles, since —E > —Vq. However, it must be recalled 
that this latter interpretation also implies that particle/antiparticle number is not conserved (only 
total charge and other additive quantum numbers are conserved) . Our spinor equation, as occurs for 
the Dirac equation, would thus no longer be a single particle equation. Some have claimed this to be 
another positive feature of the Dirac equation since it anticipated pair creation in field theory. 



III. WEYL DIUFFUSION WITH A BARRIER POTENTIAL 

The barrier potential V{x) = {Vq for < x < i , elsewhere} can be treated for calculational 
purposes as a two "step" interaction. As we have shown in previous publications [25H27] it is also 
physically correct to consider it as such. A brief argument for this is as follows. Consider a numerical 
simulation of an incoming wave packet from the left, say with E > Vo + \Jp\ + p^c so that diffusion 
occurs. Upon reaching a; = (the barrier) the wave packet splits into two amplitudes. One the 
reflected wave, the other the transmitted wave moving in the region of the potential. The group 
velocity of the latter is of course lower than the incoming/reflected wave group velocity. When the 
wave packet reaches the end of the barrier, reflection and transmission occurs anew. The very existence 
of a reflected wave for a downward step is characteristic of quantum mechanics. This procedure 
repeats itself ad infinitum. Infinite refiected and transmitted wave packets are produced. If L is large 



compared to the wave packet size, the individual wave packets will be separated and interference 
effects are negligible. Incoherence reigns. Even when the barrier size is of order or smaller than the 
wave packets this approach can be used except that now the overlapping wave packets amplitudes 
are coherent and must be summed. Separation into individual contributions may even be done for 
single plane waves which, of course, are totally coherent. The sum of the infinite contributions for 
plane waves yields exactly the same result as that deduced from a matrix calculation of the continuity 
equations. This argument can therefore be reversed. A wave packet convolution of a (matrix) barrier 
solution will yield infinite outgoing wave packets when L grows larger than the wave packet size. This 
is by no means obvious, but numerical calculations confirm this claim |25[|26) . 

We have called the above method of calculation the two step approach to the barrier. However, 
it actually requires three step results since the discontinuity at a; = receives waves impinging both 
from the left (the original wave) and from the right (reflected waves at a; = L). These three results 
arc listed below. Any two can be derived from the third by implementing the momentum changes 
indicated, together with the appropriate plane wave phase changes. 

• Step 1 [x = 0] (see section II for derivation), 

incoming momentum {E/c,Pi,p2,P3), reflected (i?/c, —pi,p2, Pa); transmitted [{E — Vo)/c,qi,p2,P3], 
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Tn = 
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Step 2 [x = L], 

incoming momentum [{E~Vq)/ c, Qi, P2,P3\, reflected [{E-Vo)/c, -(71,^2,^3], transmitted {E/c,pi,p2,P3), 



1-/3 
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• Step 3 [x = 0], 

incoming momentum [{E—Vo)/c, —qi,p2,p3\, reflected [{E—V„)/c,qi,p2TP3\, transmitted (£'/c, — ^1,^3,^3), 

f = 

^ + ^' -p,{E-V„+p3c)-tp2Vo . 

with 7 = — — ^ = (3 . 10) 

to = ' , 
1 + 7 



Observing that 1 + (3 = (1 + Q;)/Rc[a] and 1 ~ (3 — {a ~ 1)/Re[a], we can rewrite all the previous 
amplitudes in terms of a. 
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The total transmitted wave beyond the barrier {x > L) is then given by 

s t Q-iPiL/h 
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which incidently is normahzation independent. After simple algebraic manipulations, we find 
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Similarly, for the total reflected amplitude, we have 
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which results in 
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For diffusion, the single plane wave reflection and transmission probabilities are | r \ and 1 1 \ and of 
course satisfy 

\rf+ \tf = I . (14) 

The above results coincide with those calculated by solving the coupled continuity equations. However, 
it must be recalled, because often forgotten that even for diffusion \t\ is the transmission probability 
(i.e. has a physical meaning) only for a single plane wave. For wave packets small compared to the 
barrier width incoherence dominates. The transmission probability is then an infinite sum of squares 
and not the square of an infinite sum. In the limit of total incoherence the probabilities become 



Incoherent T-Probability 



Incoherent R-Probability 



IWlI 



s=0 
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(15) 
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with again probability conservation. 

Now the apparent conundrum, the step potential has only a single reflected and transmitted wave. 
How can this be reconciled with the infinite contributions (for large L)described above? The answer 
is very simple. The first contributions at x = are, obviously, the step results. The others require 
a time for their appearance which are multiples of L/vg, where Vg is the group velocity of the wave 
packet in the barrier region. As L — >■ oo this time interval goes to infinity. Formally, even for the 
step the other contributions exist but one must attend an infinite time interval to see even the second 
contribution and so forth. This feature is "hidden" when using (unconvoluted) plane waves since no 
barrier size exceeds the size of a plane wave. 

An interesting observation is that the standard resonance condition, qiL = mrh, is only valid in the 
limit of total coherence. However, another resonance condition is brought to light in the above results. 
The total transmission amplitude equals unity when there is "head on" diffusion, P2 = Pa = 0. This 
resonance condition does not depend upon the incoming energy value. Furthermore, it is independent 
of coherence requirements. The result occurs step by step for head on collision. This latter resonance 
condition is a consequence of zero mass. It does not appear for a massive Dirac particle. 



IV. KLEIN PAIR PRODUCTION 



In the previous section, we used the two step method for the calculation of the total transmission and 
reflection amplitudes. Both can be expressed as inflnite sums of amplitudes. Each term in the sum 
represents a wave packet source (once convolution is performed). The series are convergent. The total 
amplitudes can be calculated either by the two step method or by the matrix method based upon 
continuity at both x ~ and x — L. Which one uses is a matter of taste. 



Note that there is always oscillatory behavior for the y and z axis. It is only the x-axis (that of 
potential diseontinuity) where qi may be either real (oscillatory behavior) or imaginary (evanescent 
behavior). Bypassing, as previously anticipated, the evanescent case (tunneling) which exists even 
for mass zero particles when there is a transverse momentum component, we now pass to the Klein 
energy region, E kVq — ^/pfl-pl c. 

We have argued for the step result \ro\ > 1. This means that the first contribution to the total 
reflection amplitude for a barrier exceeds unity. This is already incompatible with the matrix solution 
which always satisfies |r| < 1 

The situation is even more radical. The holes (antiparticles) produced at a; = have energy —E. 
They travel within a well potential V{x) = —Vg for < x < L and zero beyond. At x = L they 
are within their own Klein zone. Thus, they also create Klein pairs, holes/electrons. The holes are 
reflected to the left. Consequently, all holes are entrapped within the barrier/well region. Pair creation 
occurs at each potential discontinuity and hence with time the hole density within the barrier/well 
increases. 

The series for t and r contain the loop factor 

_ 1+ \af -2Rc[a] .rg^L/H 

^ 1+ +2ReH 

The energy zone for diffusion, 

E>Vo + ^pl +plc =^ E>Vo-psC =^ Rc[a] > |rj^fo| < 1 • 

In this case, the series for t and r can be summed, see Ea. ([T^ and (|13p . However, in the Klein energy 
zone, 

E <Vo~ ^pl +plc E <Vo-psC =^ ReH < |rj^fo| > 1 • 

Thus, the two step series in the Klein energy zone cannot be summed. The matrix method is here 
incompatible with the two step calculation. Holes entrapped within the potential well (barrier) will 
bounce back and forth an infinitum increasing at each reflection their number. Thus, the number 
of pairs created grows without limit in time. This is of course only theoretical since we expect the 
corresponding growth in (antiparticle) charge to eventually modify the potential itself. 



IV. CONCLUSIONS 



In the previous sections, we have presented the results of the three dimensional interaction of a chiral 
(Weyl) fermion with first a step and consequently (using the step results) with a barrier. We have 
considered separately the case of diffusion, E > Vo + y/pf+p^c, and Klein zone, E <Va~ ^p\ + p\ c. 
We have bypassed tunneling in this work. 

Now for graphcne, we must set = so that our solutions apply to a physical plane (a;, y). We 
must also allow for the Fermi velocity by the substitution, c — S- • In this limit 
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It is convenient to introduce the angles 
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Thus, from Eg. dT^ we find the formulas standard to graphcne litterature [51 [T5]. 
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'sign[i?] sign[£' — Vg] — sm(j) sin 9 
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This equation clearly shows that beyond the standard resonances, qiL = mrh, there are additional 
head-on resonances for chiral fermions which are absent for massive particles. 



We now recall some of our results. 

• For diffusion, in the limit of total coherence (wave packets large if compared to the barrier width) 
the transmission probability is given by Ea. (jl2p . For wave packets small compared to the barrier 
width, incoherence dominates and the transmission probability is given by Eq. (|15p . The first term in 
the infinite sum reproduces the step result. 

• As for the Klein zone, we have reiterated our belief in Klein pair production. Obligatory, and 
generally accepted for a step potential, but not always the result quoted for a barrier potential. Our 
claim is that this is the consequence of unwittingly summing a divergent series. The coherent diffusion 
amplitudes cannot be extended to the Klein energy zone. 

• The holes created via pair production are trapped within the potential region and these localized 
holes have a continuous energy spectrum. 

The question of Klein pair production is one more hypothesis that graphene could help to test. 
If valid it implies a very different physics for "above potential" (Dirac/Weyl diffusion) and "below 
potential" (Klein diffusion) phenomena. As always we leave final judgement to experiment. 
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